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Entropy creation rate is introduced for a system interacting with thermostats (i.e. for a sys- 
tem subject to internal conservative forces interacting with "external" thermostats via conservative 
forces) and a fluctuation theorem for it is proved. As an application a time scale is introduced, 
to be interpreted as the time over which irreversibility becomes manifest in a process leading from 
an initial to a final stationary state of a mechanical system in a general nonequilibrium context. 
The time scale is evaluated in a few examples, including the classical Joule- Thompson process (gas 
expansion in a vacuum). 
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Consider "processes" between stationary out-of- 
equilibrium states of a system of particles subject 
to external nonconservative forces whose work 
is controlled by the concurrent action of ther- 
mostats. Relevant for the phenomena are objects 
like temperature, thermostats and entropy cre- 
ation. Here a "thermostat" will simply consist of 
collections of particles external to the system kept 
at constant total kinetic energy (identified with 
the thermostat temperature) by a phenomeno- 
logical force (defined by Gauss' least constraint 
principle): this is quite different from what is 
normally called a "isokinetic thermostat", keep- 
ing constant the total kinetic energy, rather than 
that of the thermostats. "Entropy creation rate" , 
(EC), is sometimes identified with "phase space 
contraction rate", (CR), i.e. the divergence of 
the equations of motion; a serious drawback is 
that it refers to the entire set of particles includ- 
ing the ones in the thermostats. However un- 
der general chaoticity assumptions a general law, 
"Fluctuation Theorem" (FT), has been derived 
governing asymptotic fluctuations of (CR): hence 
a connection between (CR) and more satisfac- 
tory definitions of (EC) is desirable. (EC) will 
be, differenly, identified with the rate of entropy 
increase of the thermostats, regarded as equilib- 
rium systems with fixed temperature as proposed 
in [1]. It is then shown to satisfy the very same 
asymptotic symmetry property expected to hold 
for the total (CR), i.e. the (FT), with the ex- 
tra feature that (FT) becomes observable on a 
time scale substantially shorter than the one nec- 
essary for observing it for all particles (for which 
it had been so far proved) including those in the 
thermostats. With (EC) and thermostats tem- 
perature it becomes possible to study processes. 
A time scale marking the waiting time necessary 
to realize the irreversible nature of a process is 
introduced here. Classical quasi static processes 
through equilibrium states turn out to have an 
irreversibility time scale +oo but the time scale 
makes sense also for non quasi static processes 



between stationary states and is discussed here. 

1. Forces and thermostats 

In studying stationary states in nonequilibrium statis- 
tical mechanics, [2, 3], it is common to consider systems 
of particles in a (finite) container Co forced by non con- 
servative forces whose work is controlled by thermostats 
consisting of other particles moving outside Co in contain- 
ers C a and interacting with the particles of Co through in- 
teractions across the walls of Co, [4]. The positions of the 
N = Nq particles in Cq and of the _/V a ones in C a will be 

denoted X a , a = 0, . . . , n, and X d = (Xo, Xi, . . . , X„). 
Interactions will be described by a potential energy 



w(x) = J2u a (x a ) + J2w a (x ,x a ) 



(1) 



i.e. thermostats particles only interact indirectly, via the 
system. All masses will be m = 1, for simplicity. 

The particles in Co will also be subject to external, 
possibly nonconservative, forces F(Xo,#) depending on 
a few strength parameters ^ = (tpi,(p2, ■ ■ ■ )■ It is con- 
venient to imagine that the force due to the confining 
potential determining the region Co is included in F so 
that one of the parameters is the volume V = \Cq\. 

Furthermore the thermostats particles will be also sub- 
ject to forces # with the property that the motions take 
place at constant total kinetic energy K a that (if ks is 
Boltzmann's constant) will be written as 



Na 1 Q 



def 3 



(2) 



and the parameters T a will define the temperatures of the 
thermostats. The exact form of the forces that have to be 
added in order to insure constancy of the kinetic energies 
should not really matter, within wide limits. 

A choice of the thermostatting forces that has been 
employed in numerical simulations has often been to de- 
fine them according to Gauss' principle of least effort for 
the constraints K a = const, [5]: the simple application 
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of the principle yields a force on the particles of the a-th 
thermostat, a = 1, 2, . . ., to be, see Eq. (1),(2), 

„q L a — U a ■ def y , . 

where i a = — 9x a W a (X , X a ) • X a is the work done per 
unit time by the forces that the particles in Co exert on 
the particles in C a . Independently of Gauss' principle it 
is immediate to check that Eq.(3) implies that the motion 
of the particles in C a conserves the kinetic energy K a . 




Caption Fig.l: The reservoirs occupy finite regions out- 
side Co, e.g. sectors C a C R 3 , a = 1,2.... Their par- 
ticles are constrained to have a total kinetic energy K a 
constant, by suitable forces # a , so that the reservoirs 
"temperatures" T a , see Eq.(2), are well defined. 

Forces and potentials will be supposed smooth, i.e. 
analytic, in their variables aside from impulsive elastic 
forces describing shocks, allowed here to model shocks 
with the containers walls and possible shocks between 
hard core particles. In conclusion the equations of motion 
of the system in its phase space T , [6] , 

T d ^ (C x R" No ) x I] (C° x B a ) EE Jo x [] Fa (4) 

a>l a 

where B a is the sphere or radius (3iV a fcsT a ) 2 in velocity 
space can be written as 

X.., i> x U • F , . j 1 Y, 

X,,, = >>x W - a j, j = l,...,N a (5) 

with i9 a given by Eq.(3) and F the external l forces in- 
troduced after Eq.(l). Denoting S*t(X, X) the configu- 
re/ 

ration into which initial data (X, X) = (X a , X a ) a= o....,n 
evolve in time t, a further condition has to be added that 
expresses the physical property that thermostats are ef- 
ficient. Namely for all initial data x — (X, X), except 
possibly for a set of data x of zero volume in T, the limit 

lim 1 f T f(S t (X,X))dt= f f(yUdy) (6) 

exists for all smooth functions / and is independent of 
(X, X), thus defining a probability distribution \x on T . 



This of course puts restrictions on the kind of forces 
acting on particles: it has to be stressed that the condi- 
tion that thermostat forces be "efficient" enough (imped- 
ing, for instance, an indefinite build up of the energy) is 
strong but it has to be assumed. It imposes on the po- 
tentials conditions that are not well understood, although 
they seem empirically {i.e. in simulations) verified with 
the simplest choices of molecular potentials, [7]. 

The notion of thermostat just described has evolved in 
the last decades and, while it will continue to evolve, it 
seems to me that it has already reached a point in which 
one can lay down the above precise notion on which to 
base a few general considerations, [6]. The class of ther- 
mostats just considered is general enough for our pur- 
poses and is amply used in numerical simulations. 

The probability distribution /z in Eq. (6) is called the 
SRB distribution or SRB statistics. Typically one is in- 
terested in studying time averages of observables depend- 
ing only on the system coordinates (X , X ) and not on 
the thermostats coordinates: in such case one needs only 
to know the distribution fi c ° obtained by restricting the 
distribution /i to such observables: also fj, c " will be nat- 
urally called the SRB distribution for the system. 

Remarks: (1) The time evolution maps St, acting on 
the full phase space T , will have the group property St • 
St> = St+t' and the SRB distribution fi will be invariant 
under time evolution. The SRB distribution \i is said 
to describe, and is identified with, a stationary state of 
the system; it depends on the parameters on which the 
forces acting on the system depend, e.g. \Cq\ (volume) 
and more generally ^ (strength of the forcings), on the 
parameters that characterize the thermostats like their 

phase space surfaces T a , i.e. (3 a d — (fcsTa) -1 , etc. The 
collection of SRB distributions obtained by letting the 
parameters vary defines a nonequilibrium ensemble. 

(2) Finite thermostats are an idealization of the appar- 
ently more appealing infinite thermostats in which the 
particles have a density p a and an empirical distribution 
which asymptotically (at infinity) is Gibbsian with tem- 
peratures T a , [6]: and in such infinite thermostats the 
temperature should stay necessarily constant because the 
local interaction between the particles in Co and those in 
the thermostats cannot alter the velocity distribution of 
the thermostat particles at large distance, at least not if 
the thermostats act as physically expected. 
(3) One could imagine replacing the forces, Eq. (3), that 
keep the temperature constant, with other thermostats 
with which the particles in the containers C a are in turn 
in contact: this is a matter of where we decide to stop in 
our investigation of the properties of what we call "sys- 
tem". Should we be interested in what really happens 
in the containers C a we would have to add them to what 
we call "system" and model their interactions with other 
thermostats. In this respect the Gaussian forces that act 
on the particles in the containers C a are a model of fur- 
ther thermostats that act on the particles in C a . If one 
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wants to avoid completely the forces (Gaussian as in Eq. 
(3) or of some other form) that fix the temperatures of 
the thermostats then one would be forced to consider in- 
finitely extended C 's. But unless the infinitely extended 
systems consist of non interacting particles, an important 
exception, see [8], not only this would lead to several con- 
ceptual problems (like existence of solutions to the equa- 
tions of motion, phase transitions in the thermostats, na- 
ture of the statistical properties of the motions inside the 
thermostats) but it would also be possibly not needed, as 
in the end we are interested only in the statistical prop- 
erties of the motions in the region Co, i.e. on the finite 
system we started with. In fact the idea of using finite 
thermostats to study nonequilibrium statistical mechan- 
ics, [5], has been the source of the progress in the field in 
the last decades. 

Remark: Besides the mentioned relation with the "ex- 
actly soluble" case of [8], the above setting is very close 
to [1] which considers what has become known as tran- 
sient properties of the fluctuations (see also [9] ) : and at- 
tempts to found on them a study of stationary states 
of systems like the ones considered here. The (CH) in 
this reference is present in the form of assumptions on 
the correlation functions on the stationary states of the 
system and thermostats. The latter works focus on ther- 
mostats which keep their temperature constant because 
they are infinite (sub)systems rather than finite isoki- 
netic (sub)systems. The finite thermostats, considered 
here, allow us to look at the problem from a different 
viewpoint. It becomes possible to define entropy cre- 
ation rate, to identify it with a contribution to the phase 
space contraction and to study the mechanism behind the 
irrelevance of the variations of the very large amount of 
energy stored in the reservoirs for what concerns the anal- 
ysis of the large fluctuation relations. In the approach of 
the present paper the condition of strict positivity of the 
average phase space contraction plays a key role and the 
restrictions on the size of the fluctuations which satisfy 
the (FT) also appears in the results; and the theory of 
the larger fluctuations (and their exponential tails) could 
be borrowed from [10]. 

2. Processes 

A process, denoted T, transforming an initial station- 
ary state fj, ini = fi under forcing <J> ini = # into a final 
stationary state /i fi n = /ioc under forcing # fi n = will 
be defined by a piecewise smooth function t — > #(i), t G 
[0, +oo) , varying between #(0) = # to ^(+oo) = #oo. 

For intermediate times < t < oo the time evolution 
x = (X, X) — > x(t) = So,tX is generated by the equations 
Eq. (5) with initial state in T and <fr(t) replacing it 
is a non autonomous equation. 

The time dependence of <fr(t) could for instance be due 
to a motion of the container walls which changes the vol- 
ume V t = \Co\: hence the points x = (X,X) evolve at 
time t in a space T(t) which also may depend on t. How- 
ever here no time dependence of the thermostats tem- 
perature will be considered. It could be introduced as in 



[1] by imagining that the potential coupling the system 
and the thermostats is time dependent so that interac- 
tion with each termostat could be switched on or off at 
will: such extra freedom and generality might require 
extensions of what is done here. 

During the process the initial state evolves into a state 
\i t attributing to an observable F t (x) defined on F(t) an 
average value given by 

(F t )= [ » t (dx)F t (x) = f [ » (dx)F t (S 0}t x) (7) 
Jr(t) Jf(o) 

We shall also consider the probability distribution /isrbj 
which is defined as the SRB distribution of the dynamical 
system obtained by "freezing" Vt, #(i) at the value that 
they take at time t and imagining the time to further 
evolve until the stationary state /J,sRB,t is reached: in 
general /i t ^ fJ>SRB.t- Of course the existence of nsRB,t 
has to be assumed a priori and this can be implied by 
the general Chaotic Hypothesis, [11]: which states that 

Motions of a chaotic system, developing on its attract- 
ing set, can be assumed to follow an evolution with the 
properties of an Anosov system. 

This means that in the physical problems just posed, 
Eq. (5) with F, i9 time independent, the motions are so 
chaotic that the attracting sets on which their long time 
motion takes place can be regarded as smooth surfaces 
on which motion is highly unstable. More precisely 

(i) around every point of T three planes M s {x), M u {x), 
M m {x) can be identified which vary continuously with x, 
which are covariant (i.e. the planes at a point x are are 
mapped, by the evolution flow St, into the corresponding 
planes around S t x) and 

(ii) the planes, called stable, unstable and marginal, with 
respective positive dimensions d s > 0, d u > and d m = 1 
adding up to the dimension of T, have the property that 
infinitesimal lengths on the stable plane and on the un- 
stable plane of any point contract at exponential rate 
as time proceeds towards the future or towards the past, 
respectively. The length along the marginal direction nei- 
ther contracts nor expands (i.e. it varies relative to the 
initial value staying bounded away from and oo): its 
tangent vector is parallel to the velocity in phase space. 
In cases in which time evolution is discrete, and deter- 
mined by a map S, the marginal direction is missing. 

(iii) over a time t, positive for the stable plane and neg- 
ative for the unstable plane, the lengths contraction is 
exponential, i.e. lengths contract by a factor uniformly 
bounded by Ce~ K l*l with C, k > 0. 

(iv) there is a dense trajectory. 

The Ergodic Hypothesis provides us with an expres- 
sion for the equilibrium averages (as integrals over the 
normalized Liouvillc distribution on the energy surface): 
likewise the Chaotic Hypothesis provides us with the ex- 
istence and a formal expression for the averages (i.e. for 
the SRB distribution), [12]. 
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Of course the hypothesis is only a "beginning" and one 
has to learn how to extract information from it, as it was 
the case with the Liouville distribution once the Ergodic 
Hypothesis guaranteed it as an appropriate distribution 
for the study of the statistics of motions in equilibrium, 
[10]. 

3. Heat generation 

The work L a in Eq. (3) will be interpreted as heat 
Q a ceded, per unit time, by the particles in Co to the a- 
th thermostat (because the "temperature" of C a remains 
constant). The entropy creation rate due to heat ex- 
changes between the system and the thermostats can, 
therefore, be naturally defined as 



a (X,X)l!/^ 



a=l 



Qa 

k B T a 



(8) 



Remarks: (1) It is natural to define Eq. (8) as "heat- 
exchanges entropy creation" rate because it is the amount 
of work per unit time that the system performs on the 
thermostats which nevertheless keep their temperature 
(identified with the kinetic energy) constant. 

(2) The above definition of "heat exchanges entropy cre- 
ation" does not require any definition of entropy itself for 
the system: it only requires the notion of entropy varia- 
tion of a reservoir considered as a system able to absorb 
work without changing it into mechanical work nor into 
its temperature variations. The latter is the "classical" 
definition of entropy variation of a reservoir. 

(3) If the thermostats temperatures are equal, i.e. K a ~ 
!-/V a /3 -1 and # = 0, or ^ is conservative so that its 
potential can be imagined merged in £/ , it is possible to 
check that the SRB distribution, if existing, is necessarily 



def 



divcrgcnce(x) 



3JV tot 

E 



ditBiix) (10) 



which has the meaning of rate at which the phase volume 
around x contracts because of the forces action. 

Phase space volume can also change because new re- 
gions become accessible (or inaccessible) so that the total 
phase space contraction rate, denoted a to t,t, in general 
will be different from erf. 

It is reasonable to suppose, and in some cases it can 
even be proved (e.g. in the example considered in remark 
(3) of Sect. 5), that at every time t the configuration So,t% 
is a "typical" configuration of the "frozen" system if the 
initial x was typical for the initial distribution hq: i.e. it 
will be a point in T(t) — V t N x R N x Y\ T a whose statistics 
under forces imagined frozen at their values at time t 
will be nsRB,t, see comments following Eq.(7). Since 
we consider as accessible the phase space occupied by 
the attractor of a typical phase space point, the volume 
variation contributes an extra of (a;) to the phase space 
variation via the rate at which the volume \F t \ contracts, 
namely 



1 d\T t \ 
I.Ftl dt 



-N 



(11) 



which does not depend on x as it is a property of the 
phase space available to (any typical) x. A calculation of 
the divergence — a T of Eq. (5), see also Eq. (3), shows 
that (keeping properly into account that the momentum 
part of phase space is a sphere) 



a r (X,X)= Y. 



3JV a -l Qa-Ua 



a>0 



3N a 



k B T a 



E 

o>0 



k B T a 



-U (12) 



^SRB 



(dx)= Z-'e-^ K ° +U ° + ^>o u « + ^«>o 



W a ) 



n a>0 s(K a 



3Na_ 

2/3 



±)-dVdX (9) 



where dx = rfVrfX = phase space volume element, 
K a are the kinetic energies of the various subsystems, 
a = 0, 1, . . . and Z is a normalization factor. Hence this 
is a Gibbs distribution in one of the equivalent, if unusual, 
equilibrium ensembles at temperature (ksft)" 1 and den- 
sities p a — j^j, a = 0,...,n. This is an important 
consistency check, following (and extending) [5] . 

4, Irreversibility degree 

It is one of the basic tenets in Thermodynamics that all 
(nontrivial) processes between equilibrium states are "ir- 
reversible" : only idealized (strictly speaking nonexistent) 
"quasi static" processes can be reversible. The question 
that is addressed in the following is whether irreversibil- 
ity can be made a quantitative notion at least in models 
based on microscopic evolution, Eq. (5). 

Write Eq. (5) as X = S(X, X) = S(x) with x = (x 1} 
. . . , ±i, . . . , x 3Ntot ), N tot = J2a N a- Given a process T let 



if U^Za 



(3N a 



~3~N — k B T ' m tne l ast s ^ e P ®(N a 1 ) has 
been neglected (for simplicity and for consistence with 
macroscopic thermodynamics where iV a are extremely 
large: it would not be difficult to keep track of it, how- 
ever). Therefore the total phase space contraction per 
unit time can be expressed as, see Eq. (3),Eq. (8), 



cr tot (X,X) = ^2 



k B T a 



N i~ u 



(13) 



i.e. there is a simple and direct relation between phase 
space contraction and entropy creation rate, [6]. Eq. (12) 
shows that their difference is a "total time derivative" . 

In studying stationary states NV t /V t = 0, and the 
interpretation of U is of "reversible" heat exchange be- 
tween system and thermostats. In this case, in some 
respects, the difference U can be ignored. For instance 
in the study of the fluctuations of the average over time 
of entropy creation rate in a stationary state the term 
U gives no contribution, or it affects only the very large 
fluctuations, [10, 13] if the containers C a are very large (or 
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if the forces between their particles can be unbounded, 
which we are excluding here for simplicity, [10]). 

Even in the case of processes the quantity U has to 
oscillate in time with average on any interval of time 
(t, oo) if the system starts and ends in a stationary state. 

For the above reasons we define the entropy creation 
rate in a process to be Eq. (13) without the U term: 

^.x)=E^-4 (») 

a 

It is interesting, and necessary, to remark that in a sta- 
tionary state the time averages of e, denoted e+, and of 
y^„ denoted V„ 9 a + , coincide because NV t /V t = 

0, as Vf = const, and U has zero time average being a 
total derivative. On the other hand under very general 
assumptions, much weaker than the Chaotic Hypothesis, 
the time average of the phase space contraction rate is 
> 0, [14, 15], so that in a stationary state: 

-' +s ?l^ i0 <15) 

which is a consistency property that has to be required 
for any proposal of definition of entropy creation rate. 

Remarks: (1) In stationary states the above models are 
a realization of Carnot machines: the machine being the 
system in Co on which external forces F work leaving 
the system in the same stationary state (a special "cy- 
cle" ) but achieving a transfer of heat between the various 
thermostats (in agreement with the second law, see Eq. 
(8),Eq. (15), only if e+ > 0). 

(2) In a stationary state atot = J2 a TTT ^ wm s& ^~ 

isfy, under the Chaotic Hypothesis of Sec.2, a fluctua- 
tion relation, see [11, 16], since the dynamics Eq. (5) 
is time reversible: the large fluctuations of the quantity 
P = i^^/o C a >o feff^ + U)dt will be controlled by 
a rate function ((p) satisfying £(— p) = ((p) — pe+, sec 
[10, 16]. However the term U will not contribute to the 
relation because U is bounded, [10]. Therefore the fluc- 
tuations of p = ji-y J Q T J2 a >o ksT dt will have exactly 
the same rate function. 

(3) Note, however, that U can be very large (being of the 
order of J2 a >o ^») ^ the thermostats are large and its 
contribution to the averages of a to t tends to as slowly 
as T~ 1 O(J2 a >0 Na)- Hence the time necessary to see 
the fluctuation relation satisfied with prefixed accuracy 
when p is defined as the average of a to t may be enor- 
mously larger than the time necessary to see the fluctu- 
ation relation satisfied with p defined as the average of 

J2 a k ®?p ■ In the case of infinite reservoirs, considered in 
[1], this leads to an "apparent" violation of the (FT) as 
pointed out in [13] quantitatively in a special case and 
extended in [10] (see the "exponential tails in the latter 
references) . 

(4) Therefore the Fluctuation Theorem holds for the 



physically interesting ^ a as a consequence of its 

validity for atot and it becomes visible over much shorter 
time scale, independent on the size of the thermostats. 
(5) The above derivation of the (FT) is somewhat unsat- 
isfactory because one cannot turn it into a mathemat- 
ically complete theorem even if one is willing to make 
strong assumptions (like the assumption that the sys- 
tem is "Anosov" in absence of external forces and of 
thermostats): this is because the lack of a boundedness 
constraint imposed on the kinetic energy in Co makes 
phase space unbounded; more physically it is the above 
"efficiency" assumption on the interaction between ther- 
mostats and system which is not properly understood 
mathematically but which is necessary in the proof. In 
this respect the results here suffer from some of the draw- 
backs in the analysis in [1] (see comments after eq. (33) 
and at the end of Sec. 4 of the latter reference) which also 
relies on an efficiency assumption on the thermostats. 
Adding a constraint that the total kinetic energy of the 
system in Co (as done in [17]) would trivially solve this 
problem but it would be physically quite unsatisfactory. 

Coming back to the question of defining an irreversibil- 
ity degree of a process T we distinguish between the (non 
stationary) state (i t into which the initial state /xo evolves 
in time t, under varying forces and volume, and the state 
HSRB.t obtained by "freezing" forces and volume at time 
t and letting the system settle to become stationary, see 
coments following Eq. (7). We call e t the entropy cre- 
ation rate Eq. (14) and ef rb the entropy creation rate in 
the "frozen" state nsRB.t- 

The proposal is to define irreversibility degree 1(F) or 
irreversibility time scale I(r) _1 of a process T by setting: 

I{T) = IP / ( (£t) ^ ~ W)sRB,t) dt (16) 

J 

If the Chaotic Hypothesis is assumed then the state Ht 
will evolve exponentially fast under the "frozen evolu- 
tion" to HsRB,t- Therefore the integral in Eq. (16) will 
converge for reasonable t dependences of V. 

A physical definition of "quasi static" transformation 
is a transformation that is "very slow". This can be 
translated mathematically into an evolution in which # t 
evolves like, if not exactly, as 

# t = *o + (l-e-T*)(*oo-*o). (17) 

An evolution T close to quasi static, but simpler for 
computing T(T), would proceed changing $ into $00 = 
<I>o+A by A/5 steps of size S, each of which has a time du- 
ration ts long enough so that, at the fc-th step, the evolv- 
ing system settles onto its stationary state at field & +kS. 
If the corresponding time scale can be taken = k^ 1 , inde- 
pendent of the value of the field so that ts can be defined 
by 5e~ Kts < 1, then I(T) = const 6~ 1 6 2 log S^ 1 because 
the variation of o-(fc+i)<5,+ — crk6,+ is, in general, of order S 
as a consequence of the differentiability of the SRB states 
with respect to the parameters, [18]. 
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5, Comments and examples 

Remarks: (1) A drawback of the definition proposed in 
Sec. 4 is that although (£t rb )sRB t ^ s independent on the 
metric that is used to measure volumes in phase space 
the quantity (et)„ t depends on it. Hence the irreversibil- 
ity degree considered here reflects also properties of our 
ability or method to measure (or to imagine to measure) 
distances in phase space. One can keep in mind that a 
metric independent definition can be simply obtained by 
minimizing over all possible choices of the metric: but the 
above simpler definition seems, nevertheless, preferable. 

(2) It is X(r) > 0, but I have not been able to prove 
(except in the limit of a quasi static process) that I(T) > 
0: this is a natural question to raise. 

(3) Suppose that a process takes place because of the 
variation of an acting conservative force, for instance be- 
cause a gravitational force changes as a large mass is 
brought close to the system, while no change in volume 
occurs and the thermostats have all the same temper- 
ature. Then the "frozen" SRB distribution, for all t, 
is given by Eq. (9) and (e srb )sRB t = (because the 
"frozen equations" admit a SRB distribution which has 
a density, see Eq. (9), in phase space). The isother- 
mal process thus defined has therefore I(T) > ex- 
cept, possibly, if it is isentropic, i.e. if initial and fi- 
nal thermodynamic entropies are equal. In fact the lat- 
ter are cases in which the initial, intermediate and fi- 
nal states are described by a density function pt(x) on 
phase space hence, [14], it is (£t rh )sRBt = and, [19], 
( £ t)u t — ~li I Pt(x) log p t (x)dx. Since the initial and 
final states are equilibrium states this means that if the 
initial and final thermodynamic entropies (i.e. Gibbs') 
are different it must be let),, ^Oona set of times of 
positive measure and I(T) > and we expect it to be of 
order S (see last comment in Sect. 4). On the other hand 
if the thermodynamic entropies are the same and T is al- 
most quasi static then the discussion at the end of Sec. 4 
shows that the irreversibility time scale may tend to +oo 
when S — > much faster than in the cases in which the 
entropies differ (e.g. as S 3, rather than as S). 

(4) The comment (3) shows that if a connection between 
thermodynamic entropy and irreversibility scale can be 
established at all it will neither be too direct nor too sim- 
ple because in the considered isothermal process linking 
two equilibrium states there is, in general, a non zero 
entropy variation hence T(T) > but, if performed very 
slowly, it can be close to be reversible both in the classical 
Thermodynamics and in the irreversibility scale senses. 

(5) Consider a typical irreversible process. Imagine a gas 
in an adiabatic cylinder covered by an adiabatic piston 
and imagine to move the piston. The simplest situation 
arises if gravity is neglected and the piston is suddenly 
moved at speed w. 

Unlike the cases considered so far, the absence of ther- 
mostats (adiabaticity of the cylinder) imposes an exten- 
sion of the analysis. In this case the phase space is the 



surface of constant energy in Cq° x R 3N ° rather than the 
full Cq" x R 3N ° . Therefore if the piston moves at a speed 
which is too slow (i.e. slower than the maximum veloc- 
ity of the particles on the energy surface) it will change 
the total energy of the gas in a way that depends on the 
special configuration initially chosen. 

Therefore the simplest situation arises when the piston 
is moved at speed so large that no energy is gained or lost 
by the particles because of the collisions with the moving 
wall (this is, in fact, a l case in which there arc no such 
collisions). This is an extreme idealization of the classic 
Joule-Thomson experiment. 

Let S be the section of the cylinder and H t — H n + wt 
be the distance between the moving lid and the opposite 
base. Let O = S H t be the cylinder volume. In this 
case the volume of phase space changes only because the 
boundary moves and it increases by N w S fi^ -1 per unit 
time, i.e. its rate of increase is Njj^. 

Hence (e t ) t is -Njjfc, while e s t rb = 0. If T = £ is the 
duration of the transformation ("Joule-Thomson" pro- 
cess) increasing the cylinder length by L, then 

l(T) = N-2£N2(jl : ) 2 dt^w 7I J^ (18) 

and the transformation is irreversible. The irreversibility 
time scale approaches as w — * oo, as possibly expected. 
If H = L, i.e. if the volume of the container is doubled, 
then I(T) = ^ and the irreversibility time scale of the 
process coincides with its "duration". 
(6) A different situations arises if in the context of (5) 
the piston is replaced by a sliding lid which divides the 
cylinder in two halves of height L each: one empty at 
time zero and the other containing the gas in equilibrium. 
At time the lid is lifted and a process T' takes place. 
In this case Vt = V6(t) because the volume V = SL 
becomes suddenly double. Therefore the evaluation of 
the irreversibility scale yields 

/•oo 

J(r') = N- 2 / N 2 5(tf dt = +oo (19) 
Jo 

so that the irreversibility becomes immediately manifest, 
I(L') = +oo, X(r') -1 = 0. This idealized experiment is 
rather close to the actual Joule-Thomson experiment. 

In the latter example it is customary to estimate the 
degree of irreversibility at the lift of the lid by the ther- 
modynamic equilibrium entropy variation between initial 
and final states. It would of course be interesting to have 
a general definition of entropy of a non stationary state 
(like the states p t at times (t 6 (0, oo) in the example 
just discussed) that would allow connecting the degree of 
irreversibility to the thermodynamic entropy variation in 
processes leading from an initial equilibrium state to a 
final equilibrium state, see [20]. 

The time scale introduced in Eq.(16) refers to the en- 
tire process rather than to what happens to the system 
only; it might have alternative interpretations: but it 
seems a quantity of interest in itself. 
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